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We study the dynamical structure factor of gapped one-dimensional spin systems in the 
critical phase in high magnetic field. It is shown that the presence of a "condensate" in the 
ground state in the high-field phase leads to interesting signatures in the response functions. 

§1. Introduction 

Recently, there has been an increasing interest in the physics of gapped one- 
dimensional spin systems subject to an external magnetic field which is strong enough 
to close the gap and bring the system into a new critical phase with finite magnetiza- 
tion and gapless excitations. ^ " The low-energy physics of the critical phase is usu- 
ally described by the Luttinger liquid, 9 )' 14 ) in some cases the problem can be directly 
mapped to the effective S = ^ chain. 7 )' 8 )' 10 )' 15 ) Such a description is equivalent to 
neglecting certain degrees of freedom (e.g., two of the three rung-triplets in case of 
the strongly coupled spin ladder 7 )' 8 )' 10 )). Those neglected states, however, form 
excitation branches which contribute to the response functions at higher energies, 
and this contribution is generally much easier to see experimentally than the highly 
dispersed low-energy continuum of the particle- hole ("spinon") excitations coming 
from the Luttinger liquid itself. At present, at least in two quasi-one-dimensional 
materials, Ni(C 2 H 8 N 2 ) 2 Ni(CN) 4 (known as NENC) 6 ) and Ni(C 2 H 8 N 2 ) 2 N02(C10 4 ) 
(abbreviated NDMAP) 5 )' 16 ) those high -energy branches were found to exhibit in- 
teresting behavior in electron spin resonance (ESR) experiments, with changes in 
the slope of the ESR lines as functions of the field. Quite recently, inelastic neutron 
scattering (INS) experiments on NDMAP, 17 ) and on TICUCI3, featuring a similar 
behavior, were reported. 

Motivated by those observations, we have studied the ESR response in the high- 
field phase of a large-D S = 1 chain. 20 ) It was shown that presence of the Fermi sea 
leads to renormalization of the energies of higher- lying excitations, which may result 
in nontrivial dependence of the excitation energy on the magnetic field. 

In the present paper we extend our previous work for the general case of arbitrary 
transferred momentum q (the ESR case corresponds to q = 0) and study the con- 
tributions of the high-energy excitations to the dynamical structure factor S(q,ui) 
in the hardcore boson approximation. We use the S = 1 Affleck-Kennedy-Lieb- 
Tasaki 19 ) (AKLT) chain and strongly coupled S = 1/2 ladder as specific examples, 
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but the general idea behind the calculation is in fact quite universal and applies 
to any gapped one-dimensional system in the high-field critical phase. We show 
that characteristic signatures in the response persist for nonzero q as well, and can 
therefore be observed in INS experiments. 

§2. Model examples and hardcore boson approximation 

In Ref. 20 ) we have used the hardcore boson approximation as an effective model 
to study the response of the large- D 5 = 1 chain. In this section, we show that the 
same type of approximation can be used essentially for any gapped ID system. 

2.1. Strongly coupled ladder 

A two-leg 5 = \ ladder with the rung coupling Jr much larger than the leg 
coupling Jl in a strong magnetic field H can be mapped to an effective 5 = \ 
XX Z chain by keeping only two states on each rung, a singlet \s) and a triplet |t+), 
and identifying them with the effective spin-i states. 7 )' 8 )' 10 ) Alternatively, one may 
describe the triplets \t+) as hardcore bosons. The resulting 5 = \ chain model 

Ws=i/2 = E J l{s%S* +1 + SlSl +1 + AS z S z n+1 } + (J R - \j L - H)S z n (2-1) 

has the anisotropy A = 1/2 of the easy-plane type and is in turn effectively described 
by the Luttinger liquid model. If the field is only slightly larger than the critical 
field H c , the density of triplets is low and the Luttinger parameter K tends to its 
non- interacting value 1. Thus, in the vicinity of H c the interaction between the 
bosons is irrelevant, and only the hardcore constraint is essential. 

Now, if one wants to include the neglected and |io) states, one can consider 
them as additional species of hardcore bosons. Neglecting all interactions between 
the bosons with the exception of the hardcore constraint, one arrives at the effective 
model of the type 

,fjPn,fj, 

n/i 

where fj, = +1, and —1 numbers three boson species (triplet components with 
S z = fi), t = Jl is the hopping amplitude which is equal for all types of particles, 
£n = zf! — I^H, and the zero-field energies are given by effl = Jr. This model is, 
of course, highly simplified, but nevertheless it captures the main physics which is 
hidden in the hardcore constraint. 

The connection between physical spin operators and bosonic operators for the 
ladder is given by the formulas S"^ 2 \ = ^{1 /2)(t\ sl +t a )—ie a p 1 Vpt 1 , where (a, j3, 7) G 

(x, y, z) and the operators t are connected to b in a standard way, b±\ = ^2~ l l 2 (t x ± 
it y ), b = t z . 

2.2. AKLT chain 

Mapping of the 5=1 chain to the effective hardcore boson model is some- 
what more complicated, but still possible. Consider 5 = 1 AKLT chain in external 
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magnetic field H described by the Hamiltonian 

T~t = |s„S n+ i + -(S n S n+ i) 2 — #5* j , (2-3) 

where S n denotes the spin-1 operator at the n-th site. The zero-field gap of the 
AKLT model is known to be A ~ 0.70 21 \ and we are interested in the high-field 
regime H > H c = A, when the gap closes. 

The ground state for H < H c is exactly known and can be compactly represented 
in the matrix product (MP) form 22 )' 23 -* 

#o = tr{ gi g 2 ...gN}, 9n = {l/Vs)(a + \-) n + a-\+) n -a°\0) n ), (2-4) 

where are the Pauli matrices in the spherical basis and \fi) n are the spin-1 states 
at the n-th site, ji = 0, ±1. 

Excitations of the AKLT model are well approximated by the triplet solitons 
(domain walls in the hidden string order), 21 ) which can be also cast into the form 
of a matrix product. 24 ) A soliton which sits at the n-th site and has S z = u is well 
approximated by the MP state 

|u,n) = ti{gig 2 ■ ■ ■ g n -i{gn^)9n+i ■ --9N} (2-5) 

Soliton states |/x,n) with different n are not orthogonal. However, one can introduce 
the equivalent set of states 

|u,n) = 2- 3 / 2 {3| / u,n-l) + | M ,n)} (2-6) 

which have the orthogonality property (//, n|//, n') = 5^^'5 nn '. The state (2-6) can 
be represented in the same MP form (2-5) with (g n o~^) replaced by the matrices f%, 

f + = \J\ i|+> - 7f - f° = y/l - 7^ + l-> - 

/" = v^/311-) + 1/V6(<7-|0) - a°\-)), (2-7) 

where the site index n has been omitted for the sake of clarity. 

One can extend this construction to many-particle states by replacing more than 
one of g's in (2-4) by /'s, however one cannot of course get the true full Fock space 
in this way, because we have four possible matrices (g, f^) for every site, while the 
original 5 = 1 problem has only three states per site. The price one has to pay is 
deviations from orthonormality for states with two or more particles close to each 
other. For example, one may observe that states |+, n; +, n'), describing two solitons 
with S z = +1, remain orthogonal, but their norm depends on the distance between 
the solitons, (+,n; +,n'|+,n; +,n') = 1 + (5/12)(-l/3)l n " n 'L The deviations 
vanish with increasing the distance between the solitons. Therefore the approximate 
description in terms of multisoliton states becomes worse as the density of solitons 
increases, which means it is useful only in the vicinity of the critical field. 

Restricting our Hilbert space to the states of the MP form (2-4) with some 
number of g matrices replaced with / + ones, one can identify the presence of the 
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AKLT matrix g and the S z = +1 soliton matrix / + at the certain site with the 
effective spin-down and spin-up states, respectively. The resulting effective S = \ 
chain is described by the following Hamiltonian 

w*=i/2 = E JeJs x n s x n+1 +sy n s y n+1 \+Y: E v m s z n s z n+m +j2(so+J2 v m-H)s z n 

n ' n m>0 ™ m 

(2-8) 

where J eS = 10/9, e = 50/27, and F m are very small, Vl = -0.017, V 2 = -0.047, 
K 3 = 0.013, V 4 = -0.0046, etc. 25 ) Thus, if one neglects the small interaction V m , 
then in the vicinity of H c the AKLT chain is effectively described by the XY model, 
i.e. by noninter acting hardcore bosons. 

Again, as in the case of the ladder, one can include the neglected degrees of 
freedom in the hardcore approximation, arriving at the model of the form (2-2) with 
e^ 0) = e and t = J e s/2. 

The connection between the physical spin operators and the bosonic operators 
for the AKLT chain is more complicated than in case of ladder. The reason for that 
is the fact that the actions of on the AKLT state are most naturally expressed 
in terms of the states n) which are in turn non- locally expressed through the 
orthogonalized states |/x,n). For example, the action of a physical spin operator Sf^ 
on the vacuum state \v) (with no bosons) is given by 

oo 

S» = (2c„/3)(&t ^ + 4 E (-m m )bl +m J |«>, (2-9) 

m=l 

where cq = ^ and c±\ = =Fl- 

§3. Dynamical structure factor in the hardcore approximation 

Thus, we are led to the problem of calculating the response functions of the hard- 
core model (2-2). The dynamical structure factor (DSF) S aa (q,Lo) in the spectral 
representation has the form 

S aa ( q ,u) = We- E f/ T \(f\S a ( q )\i)\ 2 6(co-E f + E i ), Z = E^" /T , (3-1) 
Z fi 

where S a (q) is the Fourier-transformed total spin operator. The DSF is related to 
the dynamical susceptibility x(?j through the fluctuation-dissipation theorem, 

^(^^ = l 1 _l-u ) /T ^ aa (Q,^. (3-2) 

The spin operator is generally a sum of linear and bilinear terms in bosonic operators, 
which correspond to different physical processes. The ground state at H > H c 
contains a "condensate" (Fermi sea) of b+\ bosons, and at low temperature the most 
important excitations are of the particle-hole type. We will thus take into account 
only those processes which involve states with at most one &o,-i particle: 
(A) creation/annihilation of a low-energy b + \ boson, 
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(B) creation/annihilation of one high-energy or bo) particle, and 

(C) transformation of a b + \ particle into 60 on e. 

The processes of the type (A) can be considered completely within the model 
of effective S = \ XY chain, in this case there is no need to take into account 
the high-energy states. For example, for the transversal DSF S 1 - = S xx + S yy , for 
q = 7r + k close to the antiferromagnetic wave number n one can use the known 
results 13 )' 26 ^' 27 ) for the dynamical susceptibility x xx (q,u) = x vv {,Qi uj) which yields 

r 2 f3/4 ^i/2 r(l- i^p) r(l- i^) 
X *> + k,u>) = A X (H) r Zo 7^ I ) ^ • (3-3) 

4^ r(f-^)r(f-^) 

Here A X (H) is the non- universal amplitude which is known numerically 12 ), and v 
is the Fermi velocity. A similar expression is available for the longitudinal suscep- 
tibility; 13 ) for the longitudinal DSF of the XY chain in case of zero temperature a 
closed exact expression is available as well 28 ), and for T / the exact longitudinal 
DSF can be calculated numerically. 29 ) Applying (3-2), one obtains in this way the 
contribution I A (q,ui) of the (A) processes into the dynamic structure factor. This 
contribution describes a low-energy "spinon" continuum. 

The processes of (B) and (C) type, however, cannot be analyzed in the language 
of S = 7; chain and require going back to the hardcore boson problem. 

3.1. (B)-type transitions 

Let us consider first the zero temperature case to understand the main features. 

For the case when not more than one high-energy particle is present, the model 
(2-2) can be solved exactly. 30 ) One deals essentially with a mobile "impurity" in the 
hardcore boson system; creation of the impurity leads to the orthogonality catastro- 
phe 31 ) and to the corresponding Fermi-edge type singularity in the response. 

In absence of the impurity, the eigenstates have a determinantal form 

Wkr-k N )= E (VVM)|det{^(^)}|, (3-4) 

where ipi( x ) = ~^L e%klX are the ^ ree pl ane wave functions, is the number of b + i 
particles (let us assume it to be even, for definiteness), Xj denote their positions, 
and L is the system length. The absolute value sign above (and throughout the 
paper) should be understood as a shorthand for having the antisymmetric sign factor 
attached to the determinant, which ensures symmetry of the wave function under 
permutations. The allowed values of momenta hi are given by 

h =7r+(2ir/L)/i, i = l,...,N (3-5) 

where the numbers l{ should be all different and half- integer (integer if N is odd). 
The energy of this state is E = Y^f=i{ e +i + 2icosA;j), and the total momentum 
P = J2i ki is zero (mod2-7r) when N is even. The ground state is defined by the Fermi 
sea configuration with Jj running from — (N — l)/2 to (N — l)/2, which corresponds 
to the momenta in the interval [Icf, 27r — /cp] with the Fermi momentum defined as 

k F = tt(1 - N/L) (3-6) 
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The excited configuration with one "impurity" boson having the momentum 
A can be also represented in the determinantal form 30 ) 



\(l^A)k[-k' N ) = -^J2 elP ' X °\ x o)k' 1 -k' N , Ixo}^-^ 



E 



nsn. 



det{(fi(yj)}\ . 



xo vi—vn 

(3-7) 

Here x$ denotes the coordinate of the impurity, and yj = Xj — xq are the coordinates 
of other particles relative to the impurity, the functions <fj(y) are given by 30 ) 

<Pj(y) = A 3 [e^ y+ ^ - (sin sin ^) E e^y +5 % 



where the phase shifts 8j are in our case of noninteracting hardcore particles inde- 
pendent of kj and are all equal to —tt/2, and Aj are the normalization factors. The 
functions (fj(y) were shown 30 ) to be in the thermodynamic limit L — > oo asymptot- 
ically equivalent to the free scattering states <Pj(y) — -±=e l ( k i y+Sj \ The energy E' 
and the total momentum P' of the excited state are given by 



N N 

E = E( e +i + 2t cos k j) + £ m + 2t cos A > P> = E k J + A • 

3=1 3=1 



(3-8) 



The allowed wave vectors kj and A are given by the same formula (3-5), but since 
the total number of particles has changed by one, they are different from those of 
the ground state: if N is even, the numbers Ij in (3-5) are half-integer for the ground 
state, but integer for the excited state (3-7). 



k'+q (b) 





Fig. 1. (a) The final-state configuration corresponding to (B)-type processes where a new higher- 
energy boson is added and thus the allowed wave vectors change because the total number of 
particles is changed by one; (b) schematic view of the "inter-band" transitions of (C) type which 
do not change the total number of particles. 

The matrix element (n, X\bj 1 (q)\g.s.), b^(q) = L" 1 / 2 J2 n Hi,n e '' qn -> which deter- 
mines the contribution to the response from the (B)-type processes, is nonzero only 
if the selection rules A = q, P' = P + q are satisfied, and is proportional to the 
determinant of the overlap matrix: 

((^^[■■■k'Jbliq^g.s.) = 5 x , q 5p, tP+q M fi , M fi = det{(^)} (3-9) 
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Due to the orthogonality catastrophe (OC), the overlap determinant is generally 
algebraically vanishing in the thermodynamic limit, |M/j| 2 oc L _/3 . The response 
is, however, nonzero and even singular because there is a macroscopic number of 
"shake-up" configurations with nearly the same energy. 

The value of the OC exponent (3 is connected to another exponent a = 1 — (3 
which determines the character of the singularity in the response, 

I B (q,u) = <6 /1 (-g)6t( g )) w cx 1/(lo - u^q))*, (3-10) 

where uj^(q) is the minimum energy difference between the ground state and the 
excited configuration. 

For q = ir the lowest energy excited 
configuration is symmetric against k = 
7r and is given by A = ir, kj = it ± ^-j, 
j = 1, . . . , N/2. If q deviates from 7r, 
A must follow q, and in order to satisfy 
the selection rules one has to create an 
additional particle-hole pair to compen- 
sate the unwanted change of momentum 
(see Fig. la). It is easy to see that for 
q slightly larger than it the lowest en- 
ergy is achieved when &2 = 2tt — kp 
and k\ = 2tt — kF — {q — vr); respec- 
tively, for q slightly smaller than 7rone 
has &2 = and k\ = kF + (vr — q). As 
q moves further from 7r, this configura- 
tion does not necessarily have the low- 
est energy. Indeed, as we have discussed 
in Ref. 20 ), there are also other possi- 
ble configurations with generally large 
number of umklapp-type of particle-hole 
pairs, whose energy may be lower, but 
their contribution to the response can be 
neglected because the corresponding OC exponent is larger than 1 for this type of 
configurations. 

Thus, the configuration described above (and its "shake-up" perturbations) will 
give the main contribution to the response for kF < q < 2tt — kF- The energy 
difference u>(q) with the ground state configuration which determines the position of 
singularity in (3-10) is given by 

ojfj,(q) = + 4t cos kF — 2t cos(/cf + \q — ir\), q e [kF,2ir — kp] ■ (3-11) 

For q outside the Fermi sea the symmetric configuration with the hole at k = ir gives 
the main contribution, 20 ) and the corresponding energy difference is 




H c H 

Fig. 2. The schematic dependence of "reso- 
nance" lines (peaks in the dynamic struc- 
ture factor, shown as solid lines) on the 
magnetic field. The dashed areas repre- 
sent continua. The processes responsible 
for the transitions are indicated near the 
corresponding lines, e.g. v — > — 1 denotes 
the (B)-type process of creating one boson 
with S z = — 1 from the vacuum, etc. 



Ufi(q) = e fl + 2tcosq + 2t(l + cosk F ), q i [k F , 2ir - k F \ ■ (3-12) 
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The OC exponent (5 is in both cases the same and equal to 1/2, so that the singularity- 
is of the square-root type. 

The quantity uj^(q), which determines the position of the peak in the response and 
is normally interpreted as the energy of the corresponding mode with S z = P , has a 
counter-intuitive dependence on the magnetic field: indeed, since the definition of k F 
is e+i + 2t cos kF = 0, it is easy to see that e.g. for q = it one has uj^q) = e M — e+i oc 
(1 — P )H instead of — P H as one would expect. The resulting picture of modes which 
should be seen e.g. in the INS experiment is schematically shown in Fig. 2. This 
effect may explain the ESR lines with the resonance energy proportional to 2H, which 
were observed in the high-field phase of NDMAP 16 ) (such lines would be normally 
interpreted as "forbidden" transitions with AS Z = 2). In recent INS experiments on 
the same compound, 17 ) the change of the slope in the field dependence of the q = ir 
gap for the S z = mode from zero for H < H c to 1 for H > H c was observed, which 
is also consistent with the above picture. One should have in mind, however, that 
this change of the slope from — fi to 1 — P is obtained in the model of noninteracting 
hardcore bosons, and it will be renormalized by eventually present interactions. 

At finite temperatures the edge singularity becomes damped. In the vicinity 
of the singularity, i.e. for uj — > uj^q), one can deduce qualitatively the behavior of 
S B (q,uj) from the formula derived by Ohtaka and Tanabe 32 ) for the edge singularity 
in the photoemission spectrum in case of the contact core hole potential 

where (3 = 1/2 is the OC exponent, uj = uj — to^(q), D = e + \ — 2t has the sense of 
the energetic depth of the Fermi sea, D = e+i + 2t is the width of the rest of the 6+i 
band, is defined as = 7 + Y^=2 ^ 2 "~ 1 ' > ~ 1 /? fc ~ 1 , 7 being the Euler constant 
and C(z) the Riemann zeta function, and, finally, C(q) is the (/-dependent factor 
which takes into account the modification of the overlap determinant in (3-9) 
due to the presence of "compensating" particle-hole excitations in the final state for 
q 6 [kF, 2ir — kp] (it is identically 1 for q ^ [kp, 2ir — fc^]): 

C(q) = F 2 (X 1 X 2 ,p lP2 ), F(AiA 2 , PlP2 ) = F(X lPl )F(X 2 p 2 ) - F{X lP2 )F{X 2Pl ) 

Xi = k F + 7r - q, X 2 = q, Pl = k F , P2 = n, (3-14) 

N[X] , , . N[p] 

w^nKn.'-s' 

v 11 m=l rn=l 

where k F = k F if q < vr and 2ir — k F if q > ir, and N[x] is defined as N[x] = 
L\k F — x\/2ir. Here A and P denote the discrete momenta of "holes" and "particles", 
and for P — > A one has F(X P ) — > 1. One can see that for q — > ir, as well as for q — > k F 
or q — > 2ir — k F , the factor C(q) — > 1, and it decays very rapidly if q moves away 
from those points inside the [kF, 2tt — k F ] interval. 



I B (q,u) 



C{q) (1 - p)e-M® 
4^T TJfi) 



2irT(D + D) 
DD 



p 
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3.2. Transitions of the (C) type 

So far we have been considering the (B)-type transitions. The dynamical struc- 
ture factor for H > H c will also have a contribution from (C)-type transitions cor- 
responding to the transformation of b + \ bosons into bo ones. Those processes do 
not change the total number of particles and thus do not disturb the allowed val- 
ues of the wave vector, so that there is no orthogonality catastrophe in this case. 
The action of the total spin operator S~ (q) is in this case proportional to that of 
R(q) = ^= J2 n b\b +l e iqn . It is easy to see that the action of R(q) on a state (3-4) 
characterized by the set of wave numbers {k\ ■ ■ ■ kjy} is given by 




R{qm 1 ... kN ) = x ^Y.4f^ p+q)X0 



N 



••& m _ifc m _l_i •■•k]\ 

m=l 



(3-15) 



where xo denotes the position of the created bo particle and the rest of notations is 
as in (3-7). The corresponding matrix element has a very simple form, 

((ti,\) k[ ... k > N _ 1 \R(qW{h---k N }) = jN/L6p,, P+q 6 x , km+q , (3-16) 

and the problem of calculating the response is thus equivalent to that for the ID 
Fermi gas, with the only difference that we have to take into account the additional 
change in energy Eo ~~ e +i which takes place in the transition (see Fig. lb). One 
can use the well-known formula for the susceptibility 33 -* to obtain the contribution 
I c (q,u>) = (R Ji (—q)R(q)} U j of (C) type processes into the response: 

I C (q,v) = 1 _ ~^7t ^ / dk{n[£ +1 (k)] -n[e (k + q)]}S(u-e (k + q) + e +1 (k)), 

(3-17) 

where e^{k) = e tl +2tcosk, andn(e) = (e £ / T +l) 1 is the Fermi distribution function. 
This contribution contains a square-root singularity, which survives even for finite 
temperature and is located at 



u = £o - e+1 + 2ty/ 2(1- cos q). (3-18) 

It is interesting to note that though I c (q,uj) generally does not contain any quasi- 
particle contribution, there is an exception: it transforms into a 5-function as q — > 0. 



§4. Summary 



The dynamical structure factor S(q, uj) of a gapped one-dimensional spin system 
in the high-field critical phase is studied in the hardcore boson approximation. It is 
shown that the presence of a "condensate" (Fermi sea) of S z = +1 triplets in the 
ground state in the high-field phase leads to interesting peculiarities in the contribu- 
tions to S(q,uj) from the high-energy excitations (uncondensed triplet components 
with S z = 0, —1). The energy of high-energy triplets gets renormalized because of 
the presence of the Fermi sea, which leads to a substantial change in the magnetic 
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field dependence of the triplet energies, typically from —HS Z to H(l — S z ). Extend- 
ing our previous study 20 ) , we show that such peculiarities in the response persist for 
any q, and can therefore be observed in INS experiments. 

Acknowledgments. — This work was partly supported by Volkswagen-Stiftung 
through Grant No. 1/75895. 



References 

T. Sakai and M. Takahashi, Phys. Rev. B 43 (1991), 13383. 

G. Chaboussant, Y. Fagot- Revurat, M.-H. Julien, M. E. Hanson, C. Berthier, M. Horvatic, 
L. P. Levy, and O. Piovesana, Phys. Rev. Lett. 80 (1998), 2713. 

G. Chaboussant, M.-H. Julien, Y. Fagot-Revurat, M. Hanson, L. P. Levy, C. Berthier, M. 

Horvatic, and O. Piovesana, Eur. Phys. J. B 6 (1998), 167. 

Z. Honda, H. Asakawa, and K. Katsumata, Phys. Rev. Lett. 81 (1998), 2566. 

Z. Honda, K. Katsumata, M. Hagiwara, and M. Tokunaga, Phys. Rev. B 60 (1999), 9272. 

M. Orendac, S. Zvyagin, A. Orendacova, M. Seiling, B. Liithi, A. Feher and M. W. Meisel, 

Phys. Rev. B 60 (1999), 4170. 

K. Totsuka, Phys. Rev. B 57 (1998), 3454. 

F. Mila, Eur. Phys. J B 6 (1998), 201. 

R. Chitra and T. Giamarchi, Phys. Rev. B 55 (1997), 5816. 
T. Giamarchi and A. M. Tsvelik, Phys. Rev. B 59, 11398 (1999). 
A. Furusaki and S.-C. Zhang, Phys. Rev. B 60, 1175 (1999). 
T. Hikihara and A. Furusaki, Phys. Rev. B 63 (2001), 134438. 

M. Bocquet, F. H. L. Essler, A. M. Tsvelik, and A. O. Gogolin, Phys. Rev. B 64 (2001), 
094425. 

R. Konik and P. Fendley, preprint cond-mat/0106037. 

A. K. Kolezhuk, H.-J. Mikeska, and U. Schollwock, Phys. Rev. B 59 (1999), 13565. 
Z. Honda, private communication. 

A. Zheludev, Z. Honda, Y. Chen, C. L. Broholm, K. Katsumata, and S. M. Shapiro, 
preprint cond-mat/0107416. 

Ch. Riiegg, N. Cavadini, A. Furrcr, K. Kramer, H. U. Giidel, P. Vorderwisch, and H. 
Mutka, preprint (to appear in Appl. Phys. A). 

I. Affleck, T. Kennedy, E. H. Lieb and H. Tasaki, Phys. Rev. Lett. 59 (1987), 799. 
Commun. Math. Phys. 115 (1988), 477. 

A. K. Kolezhuk and H.-J. Mikeska, Phys. Rev. B 65 (2002), 014413. 

G. Fath and J. Solyom, J. Phys.: Condens. Matter 5 (1993), 8983. 
N. Elstner and H.-J. Mikeska, Phys. Rev. B 50 (1994), 3907. 

U. Neugebauer and H.-J. Mikeska, Z. Phys. B 99 (1996), 151. 

M. Fannes, B. Nachtergaele and R. F. Werner, Europhys. Lett. 10 (1989), 633. 

Commun. Math. Phys. 144 (1992), 443. 

A. Klumper, A. Schadschneider and J. Zittartz, J. Phys. A 24 (1991), L955. 
Z. Phys. B 87 (1992), 281. 
Europhys. Lett. 24 (1993), 293. 

K. Totsuka and M. Suzuki, J. Phys.: Condens. Matter 7 (1995), 1639. 
M. Krohn, Diploma thesis, University of Hannover (2000) . 

H. J. Schulz and C. Bourbonnais, Phys. Rev. B 27 (1983), 5856. 
H. J. Schulz, Phys. Rev. B 34 (1986), 6372. 

V. Barzykin, Phys. Rev. B 63 (2001), 140412(R). 

G. Miiller, H. Thomas, H. Beck, and J. C. Bonner, Phys. Rev. B 24 (1981), 1429. 

O. Derzhko, T. Krokhmalskii, and J. Stolze, J. Phys. A: Math. Gen. 33 (2000), 3063. 

H. Castella and X. Zotos, Phys. Rev. B 47 (1993), 16186. 
P. W. Anderson, Phys. Rev. Lett. 18 (1967), 1049. 
Phys. Rev. 164 (1967), 352. 

K. Ohtaka and Y. Tanabe, Rev. Mod. Phys. 62 (1990), 929. 

See, e.g., S. W. Lovesey, Theory of neutron scattering from condensed matter. Vol. 1 
(Clarendon Press, Oxford 1984) p. 231. 



